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PACS 05. 40. -a - Fluctuation phenomena, random processes, noise, and Brownian motion. 

PACS 82.37.-j - Single molecule kinetics. 

PACS 87. 16. Nn - Motor proteins (myosin, kinesin dynein). 

Abstract. - Bartussek, Hanggi and Kissner studied a rocking ratchet system, in which a Brownian 
particle is subject to an asymmetric periodic potential together with an oscillating force, and 
found that the direction of the macroscopic current can be reversed by changing the parameter 
values characterizing the model [Europhys. Lett., 28 (1994) 459]. In this letter, we apply their 
ratchet theory to a rotary motor-protein, Fi-ATPase. In this work, we construct a model of 
a rocking ratchet in which Fi-ATPase rotates not as a result of ATP hydrolysis but through 
the influence of an oscillating force. We then study the motion of Fi-ATPase on the basis of 
molecular dynamics simulations of this coarse-grained protein model. Although in the absence of 
ATP, Fi-ATPase exhibits directionless Brownian motion when there exists no oscillating force, we 
observe directional motion when we do apply an oscillating force. Furthermore, we observe that 
the direction of rotation is reversed when we change the oscillation frequency. 
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Introduction. — Recently, theories of non- 
equilibrium statistical mechanics have been applied 
to many biological systems. Among such studies, the 
fluctuation-dissipation theorem has been investigated 
with the purpose of understanding the motion of myosins 
in an actin network [1] , the fluctuation theorem has been 
used to quantify the fluctuations of cell motion [2], and 
the Jarzynski equality and the Crooks theorem have 
been applied to RNA hairpins in order to measure their 
free energies [3]. Recent technical developments have 
made it possible to investigate small systems, including 
biological systems, and this has provided another way 
to test theories of non-equilibrium statistical mechanics 
in addition to the investigation of theoretical models of 
stochastic processes. 

Ratchet theories have been proposed to describe certain 
types of non-equilibrium statistical behavior [4-6], and 
they predict that together, spatial asymmetry and non- 
equilibrium effect, can cause directional motion in stochas- 
tic systems [5]. Among the various ratchet models, that 
consisting of a Brownian particle in an asymmetric peri- 
odic potential together with an oscillating force is called 
a "rocking ratchet." For such ratchets, it has been shown 
that the direction of the current can be reversed by chang- 
ing the parameter values characterizing the model [6] . Al- 
though the application of such current reversal to nano- 
sized systems is mentioned in Ref. [6], as far as we know, 



there is yet no such experiment. In this letter, we numeri- 
cally investigate the rocking ratchet theory in application 
to a nano-sized bio-molecule to which many theories of 
non-equilibrium statistical mechanics have been applied 
in recent years. 

We study a rotary motor-protein, Fi-ATPase, in which 
a 7-subunit rotates in the center of a3/33-subunits [7,8]. 
In this letter, we construct a model of a rocking ratchet 
in which Fi-ATPase rotates through the influence of an 
oscillating force, and particularly study the motion of Fi- 
ATPase in the absence of ATP. Because in this case no nu- 
cleotide is attached to the /3-subunits, such a state is called 
a nucleotide-free state. In the absence of ATP, it was ob- 
served that the 7-subunit exhibits directionless Brownian 
motion with rotational steps of ±120° [9]. This indicates 
that the rotational potential experienced by a 7-subunit 
is periodic with a period of 120°, reflecting the three- fold 
symmetry of ck3/33-subunits. Then, if this periodic poten- 
tial, which reflects the interaction between a 7-subunit and 
a3/?3-subunits, is asymmetric, we expect that a 7-subunit 
will display a directional rotation when we apply an oscil- 
lating force. 

In order to test this conjecture, we performed molec- 
ular dynamics simulations of a protein model, in which 
a protein molecule is regarded as consisting of a-carbon 
atoms, employing a coarse-grained representation of the 
amino-acid residues [10]. We find that the 7-subunit ro- 



K. Hayashi et al. 




Fig. 1: Schematics of the Fi-ATPase in a nucleotide-free 
state, where the Fi-ATPase consists of the a-subunits (light 
gray), the /3-subunits (gray) and the 7-subunit (black). 
{qe/3e, ob/?e, Qe/3e} is composed of the aE/3E-subunits of the 
crystal structure 1BMF [11]. 

tatcs counterclockwise when viewed from the membrane 
side in the case that the oscillation period is larger than 
the characteristic time scale of the system dynamics, and 
it rotates clockwise in the opposite case. 

Fi-ATPase in the presence of ATP. Fi-ATPase 
is a rotary motor-protein in which a 7-subunit rotates in 
the center of a3/?3-subunits due to ATP hydrolysis [7,8]. In 
the presence of ATP, a 7-subunit rotates counterclockwise 
when viewed from the membrane side and makes a 120° 
step as a result of the consumption of one ATP molecule. 
(In this letter, the counterclockwise direction is regarded 
as the plus direction.) Each step of 120° is divided into 
40° and 80° substeps, caused by the two chemical reactions 
involved in this process, i.e. the product release and the 
ATP binding, respectively. 

X-ray studies of the crystal structure of Fi-ATPase re- 
veal that a /3-subunit possesses several different struc- 
tures [11]. The ATP-bound /3-subunit, the ADP-bound 
/3-subunit and the nucleotide-free /3-subunit are denoted 
by /3tp, /3dp and /3e, respectively. While /3tp has a closed 
form at its C terminus, /3e has an open form. Because 
the structures of the /3-subunits change as the chemi- 
cal reactions involved in ATP hydrolysis proceed, it is 
widely believed that the coordinated push-pull motion of 
the C termini of /3-subunits causes the rotation of the 7- 
subunit [12]. Such substeps have been observed in numer- 
ical studies of coarse-grained protein models by changing 
the reference structures of the /3-subunits so as to create 
the push-pull motion [13,14]. In this letter, although we 
do not study such push-pull motion, we do study rotation 
in Fi-ATPase. This rotation is realized in the absence of 
ATP through use of an oscillating force in a coarse-grained 
protein model (see the Appendix). 

Fi-ATPase in the absence of ATP. While a 7- 
subunit always makes steps of +120° in the presence of 
ATP, both forward, +120°, and backward, -120°, steps 
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Fig. 2: (Top) An example trajectory of the angle, 8, of the 
7-subunit in a nucleotide-free state. 8 = 0° represents the 
position of the 7-subunit in the crystal structure 1BMF [11]. 
The dotted lines represent 8 = 20° + 120°n (n = 0, ±1, ±2 • • •). 
(Bottom) The rotational potential, — log(P(#)). (Inset) The 
time correlation function, (9(t)0(0)) — (8(0)) 2 , in dwell time 
(circles). The fitting curve (a thick black curve) represents 
2.2exp(-t/5,236At). 



can be observed experimentally in the absence of ATP 
[9]. Because the structure of /3e, which has an open form, 
makes the height of the rotational potential smaller [9], 
both +120° and —120° steps can be observed. 

First, we investigated the Brownian motion of Fi- 
ATPasc in equilibrium (without an oscillating force) 
through numerical simulations. Here, we assume that the 
reference structure of the a3/?3-subunits in a nucleotide- 
free state can be represented by {o;e/3e, «e/3e, cxePe}, and 
that in equilibrium there is no structural change of the /3- 
subunits. We computed the time evolution of the rotation 
angle 9(t), which represents the direction of the a-carbon 
atoms in the 7-subunit encircled in the top-left region of 
Fig. 3. (See the Appendix and Ref. [13] for a precise def- 
inition.) In the upper graphs of Fig. 2, we plot a typical 
trajectory of 9(t) in the nucleotide-free state as a function 
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Fig. 3: (Top) Schematics of the rocking ratchet based on Fi- 
ATPase in a nucleotide-free state. (Bottom) In the case that 
there exists an oscillating force with tq — 10, OOOAt and A = 
0.005, the 7-subunit exhibits counterclockwise rotation. The 
snapshots were taken using VMD [15]. 

of time. Directionless Brownian motion is observed, as 
there are both steps of +120° and —120°. In the upper 
graphs of Fig. 2, it is seen that the 7-subunit stops at 
angles of approximately 20° + 120°n (n = 0, ±1, ±2, • • •). 
We note that it has also been observed experimentally 
that in the absence of ATP, the 7-subunit stops at angles 
shifted by —20° from the ATP binding angle [9]. Because 
the ATP binding angle is considered to be 8 — 30°-40° in 
our model [13], the stopping angles plotted in the upper 
graphs of Fig. 2 seem to correspond to the experimental 
results [9]. 

In order to investigate the system dynamics in the intra- 
well of the rotational potential, we analyzed trajectories 
of 8(t) that is fluctuating around an angle of 20° + 120°n 
(n = 0, ±1, ±2 ■ • ■). By using the probability distribution, 
P(0), we compute the rotational potential — log(P(#)). 
Then, its asymmetric profile is seen in the lower graphs 
of Fig. 2. In the inset of the lower graphs of Fig. 2, 
we further plot the time correlation function, (9(t)6(0)) — 
(8(0)) 2 . From the fitting curve, the intra-well relaxation 
time is estimated to be 5, 236 At. 

Rocking ratchet. — Having verified that our protein 
model appears to be capable of describing the motion of 
Fi-ATPase in the absence of ATP, we next investigated 
the case in which an oscillating force is applied to the 7- 
subunit in order to create a rocking ratchet. This force 
mimics one that can be applied by a magnetic tweezers 
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Fig. 4: Example trajectories in the case of an oscillating force 
with the cases m = 50, OOOAt (Top) and rn = 4, OOOAt (Bot- 
tom) for A = 0.005. The dotted lines represent 6 = 20° + 120°n 
(n = 0,±l,±2--.). 

(Fig. 3). Explicitly, we apply the oscillating force 

(fx,fy) = (Asin27rt/rn sin 6*, — Asin27ri/Vn cos 6*) (1) 

to the a-carbon atoms in the 7-subunit that are not inside 
the a3/?3-subunits (marked in red in Fig. 3). Note that 
although in real experiments an oscillating force should be 
applied to the bead to which a 7-subunit is attached, for 
simplicity, here we apply this force directly to the a-carbon 
atoms in the 7-subunit. In the lower graphs of Fig. 4, we 
plot typical trajectories in the case with an oscillating force 
of periods tq = 50, OOOAt and tq = 4, OOOAt, respectively. 
In Fig. 4, it is seen that when A = 0.005, the range of 
oscillation of the 7-subunit is approximately 120° which 
is the period of the potential, and that the 7-subunit is 
likely to stop at the same angles as in equilibrium even in 
the case with the oscillating force. 

In Fig. 5, we plot 20 trajectories in the cases tq = 
50, OOOAt and m = 4, OOOAt. There, we can see clearly 
that rotation has a directional preference. Interestingly 
the average direction of rotation in the case with tq = 
50, OOOAt is opposite to that in the case tq = 4, OOOAt. It 
is thus seen that current reversal occurs when we change 
the oscillation period. Comparing these results with the 
directions of rotation found in Ref. [6], they lead us to 
conjecture that the slope of the potential in the system 
studied here, like that studied there, is somewhat smaller 
in the counterclockwise direction (the plus direction) than 
that in the clockwise direction. Indeed, the rotational po- 
tential plotted in the lower graphs of Fig. 2 is consistent 
with this conjecture. 

In order to determine the value of at which the cur- 
rent reversal occurs, in Fig. 6, we plot the average currents 
(®( T f)/ T f) as functions of tq for the several values of the 
oscillation amplitude. (Here, 17 = 2 x 10 6 At.) We find 
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Fig. 5: Typical trajectories in the cases in which the oscillating 
force has periods, tq = 50, OOOAt (Top) and r n = 4, OOOAt 
(Bottom) for A = 0.005. 

that the current vanishes at a value Tq ~ 6, OOOAt. Al- 
though we do not yet have a complete understanding of 
the mechanism of the current reversal, we believe that rj^ 
is related to the intra-well relaxation time in equilibrium, 
which is estimated to be 5, 236 At (see the inset of the 
lower graphs of Fig. 2). In Ref. [6], it is asserted that the 
intra-well relaxation time determines the oscillation pe- 
riod, which causes the current reversal in a rocking ratchet 
system. Taking these into consideration, the current re- 
versal observed in our protein model seems to be similar 
to that phenomenon originally reported in Ref. [6] . 

Discussion. — In conclusion, we investigated motion 
of Fi-ATPase in the absence of ATP by performing numer- 
ical simulations of the coarse-grained protein model, and 
fabricated a rocking ratchet applying the ratchet theory 
studied in Ref. [6] to Fi-ATPase. We observed the current 
reversal when we changed the oscillation frequency and 
predicted the directions of rotation when there exists an 
oscillating force. We hope that our predictions for a rock- 
ing ratchet consisting of Fi-ATPase in a nucleotide-free 




Fig. 6: The average currents as functions of the oscillation pe- 
riod, rn, in the cases A = 0.0025 (squares), A — 0.005 (circles) 
and A = 0.0075 (triangles). 

state will be confirmed experimentally, and that combi- 
nation of theories of non-equilibrium statistical mechanics 
with protein simulations will yield new understanding of 
statistical properties possessed by proteins [10,16]. Below, 
we discussed two things related to our results. 

It should be noted that the parameter values used in 
our simulation do not correspond to room temperature 
and the real time scale of experiments. The temperature 
used in our simulations corresponds to one third less than 
the room temperature in order for the 7-subunit not to 
lose its shape when the oscillating force is applied. How- 
ever, we believe that the basic prediction of our model - 
that the direction of rotation can be reversed by altering 
the period of the oscillating force - would be unchanged if 
such parameter values were used. In previous numerical 
studies [13, 14] also, temperatures much less than room 
temperature were used in order to make the coupling be- 
tween the structure change of the /3-subunits and the rota- 
tion of the 7-subunit "tight," on the basis of the idea of the 
push-pull model [12]. However, it is worth simulating the 
motion of Fi-ATPasc with realistic parameter values that 
represent the physiological conditions in which the molec- 
ular softness and thermal fluctuations would be more im- 
portant. With such conditions in mind, the ratchet-based 
mechanism has been considered instead of the push-pull 
model [17] in an attempt to account for the results of a 
recent experiment [18]. The same idea has been applied 
to an F -motor system [17]. 

One reason for investigating systems like that consid- 
ered here is to gain an understanding of motor proteins 
that could lead to the design and fabrication of new nano- 
machincs that can operate efficiently in thermal noise. For 
this purpose, it is important to be able to cause directional 
motion of proteins not through use of ATP hydrolysis but 
by other energy sources. For this reason, we have stud- 
ied a system in which the fluctuations and softness of Fi- 
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ATPase are exploited in such a manner to cause rotation 
through the influence of an oscillating force. Because the 
ratchet theories [4-6] are regarded as providing a theoret- 
ical foundation for the design of nano-machines, we hope 
that this study provides a basis for the development of one 
class of such machines. 

Appendix. — Here we describe the model studied in 
this work. The position and velocity of the ith a-carbon 
atom in Fi-ATPase are denoted by and Vi. The Hamil- 
tonian is given by H({r h Vi}) = J2iJT m\ v i\ 2 / 2 + V a 3 ^ 3 + 
V~t + Vint + Vfi x , where m is the mass of the ith a-carbon 
atom and N and n are the number of a-carbon atoms in 
the a3/?3-subunits and in the 7-subunit, respectively. The 
quantity V a3 p 3 is the elastic network potential for the a^(3^ 
subunits and is given by 

AT-l , 

V as p 3 ({ri} a3 p 3 ) = ~ r ' t+l1 ~ l r i ~ r ^i |)2 

i=l 

+ Ey([rrr w |-|ff-a 

N-l N 

+ E E ^(In-rA-W-rW, (2) 

i— 1 j—i+3 

where {r } represents the reference structure, and we have 
k = when \rf — r?| > r c . The elastic network poten- 
tial, V-y, is defined similarly to V a3 p 3 . Here, {f°}a 3 /3 3 is 
composed of the aE/?E-subunits of the crystal structure 
1BMF [11] , which was taken from the RSCB Protein Data- 
Bank [19]. We rotate them by ±120° in the xy-pl&ne and 
choose the reference structure as {cv^Pe, cte/3e; &e/3e}i 
where the xy-plane is defined so that the three /3g are 
on the plane. The set {t°} 7 is also constructed by using 
1BMF. The interaction potential between the 7-subunit 
and the a3/?3-subunits consists of a repulsive potential and 
an electrostatic potential: 

N N+n ( , s 12 

vwm) = E E {- 

i=l j=N+l K v 1 Jl 7 

|n-r 3 -| 2 J 

Here, % represents the sign of the electric charge of the 
amino acid around the zth a-carbon atom: qi = 1 when 
the amino acid is lysine or arginine, g,- = — 1 when it is 
glutamic acid or aspartic acid, and = otherwise. In 
our model, 7272 and the three (3g are fixed by including 
the potential Vfi x (r,) = kh\rt - rf\ 2 /2. 

The time evolution of the zth a-carbon atom is described 
by an over-damped Langcvin equation with friction coef- 
ficient r and temperature T. In our simulation, the pa- 
rameters are set to be m = 1, k = 10, N = 2823, n = 122, 
r c = 10, D = 6, T = 0.01, r c = 10, T = 0.3, e x = 3.0, 
£2 = 0.05 and kh = 0.5. We chose these values so that 



the results obtained without the oscillating force are con- 
sistent with the results observed experimentally [9]. 

Regarding the definition of the rotation angle 9(t), as is 
done in the previous study [13], we consider the projection 

X xy (t) Of X(t) = E je7 79-90^W-^78(t)]/|n(t)-r 778 (t)| 

onto the xy-planc. 9(t) is defined by the angle between 
X xy (t) and \ x ® MF where \ x ® MF is the corresponding vec- 
tor calculated for 1BMF [11]. 
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